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Abstract 

O, 

I show that a lattice theory of massive interacting fermions in 2n + 1 dimensions 

may be used to simulate the behaviour of massless chiral fermions in In dimensions if 

OS 

-t— > 
03 



the fermion mass has a step function shape in the extra dimension. The massless states 
arise as zeromodes bound to the mass defect, and all doublers can be given large gauge 
invariant masses. The manner in which the anomalies are realized is transparent: apparent 
chiral anomalies in the 2n dimensional subspace correspond to charge flow into the extra 
dimension. 
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A successful chiral regulator for gauge theories has proven elusive. Dimensional reg- 
ularization explicitly violates the chiral gauge symmetry, as does the conventional Pauli- 
Villars technique; zeta function regularization fails because the chiral Dirac operator does 
not provide an eigenvalue problem. A lattice regularization technique would be desirable, 
allowing one to simulate the standard model as well as strongly coupled chiral gauge the- 
ories [0. However, the naive approach suffers from the existence of doublers that render 
the theory vector-like instead of chiral . Some more sophisticated approaches appear to 
suffer other obstructions ||-[7|, although some have yet to be thoroughly tested [||. Any 
such regulator must correctly reproduce the anomalous Ward identities of the continuum 
theory. Taking the continuum limit must fail in theories where the gauge currents have 
anomalous divergences; on the other hand, nonconservation of global currents must be al- 
lowed in the continuum. As Banks has emphasized recently, a necessary criterion for chiral 
lattice theories is that they must correctly describe proton decay in a simulation of the 
standard model [0. The point is that the standard model has anomaly free chiral gauge 
currents, but anomalous global currents, such as baryon number 0. A lattice regulator 
that achieves this while preserving gauge invariance will have to be devious, as "no-go" 
theorems imply ||10|| . 

In this Letter I show that chiral fermions in 2n dimensions may be simulated by Dirac 
fermions in 2n + l dimensions with a space dependent mass term. In particular, when there 
is a domain wall defect in the mass parameter, the lattice theory will exhibit massless chiral 
zeromodes bound to the 2n dimensional domain wall, as well as doublers with the opposite 
chirality. Since the original 2n + 1 dimensional theory is vector-like, the doublers may be 
removed by means of a gauge invariant Wilson term. This system can be used to simulate a 
2n dimensional chiral gauge theory (for n — 1, 2) provided that the massive fermion modes 
which are not bound to the mass defect decouple, and provided that the gauge fields can 
be constrained to propagate along the defect. I argue that whether or not these conditions 
can be met depends on an algebraic relation which is just the one that ensures that the 2n 
dimensional gauge currents are exactly divergenceless. Global currents, however, may have 
anomalous divergences even in the regulated theory before taking the continuum limit, due 
to the curious role played by the extra dimension. 

I begin by recapitulating some well known facts about chiral zero modes of a free 
Dirac fermion in In + 1 dimensions coupled to a domain wall. I work in Euclidian space 
and take the 2n + l coordinates to be z M = {x, s} where a runs over 1, . . . , 2n + 1 while x is 



the coordinate in the 2n dimensional subspace. I use hermitian gamma matrices satisfying 
{Tm»7^} = 2( W and define 

7 2ri+ i = r 5 . (i) 

Note that r$ is the chiral operator in the 2n dimensional theory. The fermion mass is 
assumed to have the form of a domain wall — a monotonic function of s with the asymptotic 
form 

m(s) — ► ±m±, m± > . (2) 

s — »±oo 

The Dirac operator is then given by 

K D = [0 + m(s)] = [7 ■ V + T 5 d s + m(s)} . (3) 

Since the mass m(s) vanishes at s = one might expect massless states bound to the 
defect. Such "zeromodes" would have to satisfy Kd^q = 7 • V^o i n order to describle 
massless propagation along the defect. Specifically, the zeromodes must have the form 

V± = e^(P ± (s)u ± (4) 

where u± are constant 2 n component chiral spinors 

T 5 u± = ±u± 

and the functions (j)±(s) satisfy 

[±d s + m(s)](P±(s) = . (5) 

The solutions to the above equation are 

(p±(s) =exp ( =f / m(s')ds') , (6) 

however with the asymptotic form for m(s) given in eq. @, only <fi+(s) is normalizable. 
Thus there is a single, positive chirality massless fermion bound to the mass defect. 

What makes this model of interest for a lattice approach to chiral fermions is that the 
2n + 1 dimensional theory has neither chirality nor anomalies, yet possesses a low energy 
effective theory (well below the asymptotic mass scales m±) which is a theory of a massless 
chiral fermion in 2n dimensions. As I will discuss below, the anomalies in the effective 2n 



dimensional theory are realized through finite Feynman diagrams and therefore have their 
canonical form even in the regularized theory. 

A similar model can be formulated on a 2n + l dimensional lattice with lattice spacing 
a and sites labelled by z = n z a. The action is given by 

So = J2 ^zK D ^ z (7) 

z 

where the lattice Dirac operator is given by 

2n+l 

2a 



1 2n+l 

K D ^ Z = ^-J2 7m(**+A - **-a) + m ( s )** (8) 



i*=i 



and ft corresponds to a displacement by a in the z M direction. For simplicity I take the 
fermion mass to be a step function 



m(s) = m 9(s) = smh (°^°) fl( s ) ? d (s) = { 0, .s = . (9) 



— 1, s < —a 
0, s = 

+1, s > +a 



Then there are two chiral zeromode solutions \I/ given by 

*J(p,z) = e 4 **0±( S ,j?Kt (10) 



where the transverse wavefunctions are given by 

-jt*o|«| 



<P+{s,P) =e 



(11) 



These solutions satisfy 

K D ^(p,z) = (i/a)^-sm(pa)y+(p,z) , r 5 ^± = ±tf± (12) 

and, unlike the continuum example, are both normalizable and localized along the mass 
defect at s = 0. For \pa\ <C 1, ij ■ sin(pa)/a — > ip 1 , the inverse propagator for a massless 
mode travelling in 2n dimensions. Therefore the mode ^q for small p corresponds to the 
positive chirality state ([|) found in the continuum. However, the lattice model is seen to 
actually describe 2 2n+1 massless modes in the continuum, rather than one, corresponding 
to \1/q" and i&Q with p near the 2 2n corners of the 2n dimensional Brillouin zone. The 
conventional analysis of doubler modes @ reveals that these 2 2n+1 modes correspond to 



2 n positive chirality states and 2 n negative chirality states. The resultant lattice spectrum 
is therefore completely vectorlike in the continuum, unlike what I wish to describe. 

Nevertheless, since the 2n + 1 dimensional theory (0) is itself vectorlike, it is possible 
to add a Wilson term to eliminate the unwanted doublers. Thus I add to Sd the action 



z 
2n+l 2n- 



2n+l 



(13) 



n=i 



/*=i 



Zeromodes will now correspond to solutions satisfying eq. (|i~2|) with Kd replaced by 
Kd + awKw- Thus the transverse wavefunctions 4>± must satisfy 
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±d s + moO(s) + waA s + (2w/a) y^(cos(pid) — 1) 



i=l 



i±(a)=0 



(14) 



I will analyze this equation for the particularly convenient choice of w = .5, though there 
is actually a broad range of values that will have a similar solutions. Defining 



2n 



F(p) = ^(1 - cos apt) > 



i=i 



eq. (|TJ) becomes (for w = .5) 



±(s ± a) = -m eff (s)0±(s) , 

+amo — -F(p ) — 1 s > 



m e ff(s) = am 6(s) - 1 - F(p) 



(15) 



-amo — i^(p) — 1 s < 



In order for there to exist normalizable solutions for 4>±(s,p), |ra e fj| must satisfy the 
following conditions: 



0-1- : |m e ff | > 1 for s < 0, \m e ff\ < 1 for s > 0; 

0_ : |m e ff | < 1 for s < 0, |m e ff| > 1 for s > . 



(16) 



Since amo and F(p) are positive, the conditions for a normalizable </>_ solution can never 
be met, and so all 2n of the ^q modes have been elir 
are normalizable (f>+(s,p) solutions for all p satisfying 



be met, and so all 2n of the \& modes have been eliminated. On the other hand, there 



< am - F(p) < 2 . 
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(17) 



Thus, for < amo < 2 there is always a normalizable solution for small p, since F(0) = 0. 
Furthermore, on the Brillouin boundary F > 2, and so all of the doubler solutions for <f>+ 
have been removed for the same range of amo ■ This shows that the Wilson action ([13]) has 
eliminated all but one of the 2 2n+1 zeromodes bound to the domain wall in the continuum 
limit. Only a single positive chirality mode remains. 

It may seem surprising that the Wilson term could remove an odd number of chiral 
modes, since one expects them to be removed in pairs with opposite chirality. In this case 
though, the chiral doubler modes are paired with modes "at infinity", rather than with 
each other. The effect is more easily understood in finite volume. Consider a lattice of 
length L in the s direction, with periodic boundary conditions. Now if there is a domain 
wall at s — 0, there is necessarily an "anti-domain wall" at s — ±L/2. For every zeromode 
at s = 0, there is a zeromode of opposite chirality at s = ±L/2. The vectorlike nature of 
the 2n + l dimensional theory is now manifest: every mode is chirally paired, though the s 
dependent mass function has allowed one to separate some chiral pairs by a distance L/2 
with exponentially small overlap. 

At finite volume one can also see how the peculiar discontinuity in the spectrum 
implied by eq. ([IT]) comes about. At finite L with periodic boundary conditions, only 
the p = mode is exactly chiral; for p > satisfying the condition (|17|), the modes have 
masses exponentially small in L, corresponding to a finite probability that the mode can 
tunnel to the anti-domain wall and flip chirality. As L — > oo, this probability goes rapidly 
to zero and the modes are exactly chiral as seen above. For values of p not satisfying eq. 
( |17D the modes are massive for any L, finite or infinite. 

I conclude that at infinite volume the spectrum of the free theory contains (i) doublers 
with mass 0(w/a); (ii) Dirac fermions with mass O (mo) that propagate in the full 2n + 1 
dimensions; (iii) an exactly massless chiral zeromode bound with r§ = +1 bound to the 2n 
dimensional defect. For momentum \p\ -C 1/a the mode has the conventional continuum 
propagator. The theory is vector-like at every step: both the doublers with mass at the 
cutoff, as well as the mass mo scattering states, are all Dirac fermions. The zeromodes are 
also paired: in finite volume one sees that there are zeromodes of opposite chirality bound 
to an anti-domain wall, which are carried off to spatial infinity in the infinite volume limit. 

When background gauge fields are coupled to the 2n + 1 dimensional Dirac theory, 
there appears to be a paradox. One would expect that at energies well below mo all massive 
fermion degrees would decouple, yielding an effective 2n dimensional theory of zeromodes 
at one of the domain walls. However, such a theory is in general anomalous and cannot 
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be simultaneously gauge invariant and chirally invariant. In the continuum, the zeromode 



gauge current J zrn should have a divergence equal to [pA| - fL2 

where 

C n = (-l) n /(2 2n ir n nVj , 

F^ = ±1 is the chirality of the zeromode, and the indices run over the In dimensions. On 
the other hand, the lattice theory is manifestly finite and gauge invariant. So what is going 
on? 

The resolution of the paradox is that the massive fermion states that live off the 
domain walls do not entirely decouple — even in the low energy limit — and so the effective 
theory is not truly In dimensional. Callan and Harvey first described this effect some 
years ago |T3| while studying the properties of zeromodes on topological defects and their 
relation to the anomaly descent relations [ITJJjTH]. They showed that the combined effect 



of the background field strength F a p and the space dependent mass term 

Cm = +m(s)** (19) 



was to distort the heavy fermion vacuum, producing a Goldstone-Wilczek current [[16| far 
from the domain wall, given by 

^ 7774 S } 

{i'lfi^GW = -» if U C n^ij,a 1 ...a 2n F aia2 ■ ■ ■ -F Q2n _ lQ2 „ • (20) 

Note that the factor m(s)/\m(s)\ — 9(s) for a domain wall mass, and that this legacy from 
the heavy modes does not vanish in the limit |rao| — > oo 

Because 9(s) has opposite signs on the two sides of the domain wall, the current 
has a nonzero divergence: in the presence of nonzero E and B fields, the current has a 
net flux onto or off of the domain wall. Callan and Harvey pointed out that this flux 
precisely accounts for the anomalous divergence ([H|) of the zeromode current. Therefore, 
what appears as an anomalous divergence in the effective 2n dimensional theory is simply 
charge flowing on or off the mass defect from the extra dimension. Note that the sign of the 
mass coupling in eq. (|19|) determines both the sign of the Goldstone-Wilczek current ( |20| ) 
and the chirality T$ of the zeromode on a given mass defect, which figures in the anomaly 
equation (|I8|) . Callan and Harvey also showed that the current (|20| ) may be directly 



computed from an effective action, which is just the Chern-Simons term [I7J proportional 



to m{s)/\m{s)\ that was induced in the effective theory by integrating out the massive 
fermion modes |18| 111. 



The extra dimension is the loophole in the Nielsen-Ninomiya theorem [[10] through 
which the fermions have wriggled. The currents of the 2n + l dimensional lattice are in fact 
exactly conserved. Nevertheless, there is anomalous current violation in the In dimensional 
effective theory due to charge flowing in from the extra dimension. As the Chern-Simons 
action is the only new local operator that can be induced in the effective theory that does 
not vanish as the fermion mass increases, the only effect the massive fermions have on 
the low energy theory (aside from renormalizing couplings) is to correctly reproduce the 
anomalous Ward identities. 

As the Goldstone-Wilczek current ( p0|) does not vanish in a regulated theory, the In 
dimensional anomalous Ward identities may for the first time be directly measured on 



the lattice. Such a numerical computation has recently been performed by K. Jansen 1 20f| 



whose results show that starting from 2+1 dimensions, one can measure the 1+1 anomalous 
divergence of a chiral current. He finds that without the Wilson term ([T3p the theory is 
vectorlike and all 2n dimensional flavor currents are conserved. However, for nonzero 
Wilson coupling w, he finds a nonzero divergence for the currents, and the rate of charge 
violation on the domain wall is proportional to the applied electric field, as prescribed 
by the anomaly equation ([18]). Chiral theories with nontrivial anomaly cancellation can 
also be simulated in background gauge fields. The simplest example of such a theory is 
the 3-4-5 chiral Schwinger model: one starts with a 2 + 1 dimensional theory with three 
Dirac fermions \&3, W4 and ^5 coupled to a U(l) gauge field with charges q = 3, 4 and 5 
respectively. The mass term for these fermions is given by 

£ m = m( S )[* 3 ^3+^4^4-*5^5] , (21) 

so that the effective 1 + 1 dimensional theory on the s = defect consists of charge 
q = 3,4 fermions with F^ = +1, and a q = 5 fermion with F^ = —1. As the 2-dimensional 
anomaly in the gauge current is proportional to Trq 2 F$, it vanishes; the Chern-Simons 
term induced by the heavy fermions is proportional to Tr q 2 (m/\m\), and so it and the 
resultant Goldstone-Wilczek current also vanish. As has been confirmed numerically by 



See Naculich [IS] for a discussion of the subtleties involved. 



Jansen [|20], the gauge current is anomaly free, even though the individual flavor currents 



are anomalously divergent. 

The next step is to introduce dynamical gauge fields. Including the usual 2n + 1 
dimensional Yang-Mills action is not appropriate for describing a 2n dimensional gauge 
on the domain wall. There would exist an additional polarization for the gauge bosons, 
corresponding to the 2n + l component of the vector potential and appearing as an adjoint 
pseudoscalar in the 2n dimensional world. Another problem is that the gauge bosons 
would propagate off the defect, giving rise to the wrong force law between the zeromodes. 
The conventional 2n + l dimensional gauge action is not required though, since the action 
need only be invariant under gauge and 2n dimensional Lorentz transformations. Thus I 
consider a gauge action of the form 

2n 

S g =lh J2^U 4 \ Dtj +2(3 2 J2^U 4 \ Dts . (22) 

In general, (3\ and (3 2 are independent and can be functions of s. It appears that simply 
setting to zero the (renormalized) coupling f3 2 would suffice for simulating a 2n dimensional 
gauge theory at each s-slice. However, in the interacting theory, this is clearly not always 
possible. It is instructive to consider the continuum Yang-Mills equations for these gauge 
bosons in 2n + 1 dimensions, 

IhDjFij + D s (3 2 F ls = J z 

(23) 
faDiFgi — J s . 

From the second equation one sees that F is ~ J s /(3 2 - Therefore, taking the limit (3 2 — > 
does not eliminate the unwanted F{ 8 term from the first equation — unless J s vanishes. 

There are two possible contributions to J s , from the zeromodes and from the 
Goldstone-Wilczek current (|20|). In fact in the infinite volume limit there is no zeromode 
contribution to the to J s since / i/>7 2n+1, i/> = i/jV^ = 0, where ip(x) is any of the zeromodes 
on the mass defectH. Therefore, for J s to vanish, the component of the Goldstone-Wilczek 
current orthogonal to the domain wall must vanish. As we have seen, this occurs if, and 
only if, the gauge currents in the 2n dimensional theory are anomaly free. 



At finite volume there is an exponentially small probability for the zeromode to tunnel to 
the anti-domain wall, contributing to J s - 
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I conclude that if the 2n dimensional gauge currents are anomalous, there is nothing 
one can do to make the theory look like a In dimensional gauge theory — it looks inherently 
2n + l dimensional, where it makes sense. However, I speculate that for anomaly free chiral 
gauge currents, the fo - *■ limit decouples the unwanted gauge degrees of freedom, and 
the lattice model looks like a In dimensional chiral gauge theory. In that case, 3- and 5- 
dimensional models of this nature are expected to possess second order phase transitions, 
corresponding to continuum chiral gauge theories in 2- and 4-dimensions. Note that the 
effective In dimensional gauge coupling will be related to f3\ by a factor of the normalization 
of the transverse wavefunction + (s), correctly accounting for the change in dimensionality. 

Suppose the mechanism described here works and one can simulate the standard 
model — how does the proton decay? As I argued above, the only requirement for recovering 
a 4-dimensional gauge symmetry was that the chiral gauge current be anomaly-free. The 
global currents, such as B and L will still be anomalous. Given an SU(2) instanton event 
in the s = mass defect, a proton will flow off the domain wall and its baryon number will 
reside in the vacuum polarization of the extra dimension. Since B — L is anomaly free, 
the vacuum will disgorge a positron into our 4-dimensional world (I am assuming a single 
family). As far as we are concerned, the proton has decayed; little do we suspect that it is 
hovering nearby but out of sight in the extra dimension. 
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